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The PHENIX Run 4 Au+Au dataset provides a powerful opportunity for exploring the 
angular anisotropy of identified particle yields at high pt- Complementing traditional V2 
measurements, we present 7r° yields as a function of angle with reaction plane, up to pt 
10 GeV/c. The centrality dependence of the angular anisotropy allows us to probe the 
density and path-length dependence of the energy loss of hard-scattered partons. We will 
discuss various mechanisms for particle production in this high pt region. 

1. Introduction 

One of the most intriguing puzzles in RHIC physics is the origin of the azimuthal 
anisotropy of particle yields at high pt (> 5 GeV/c) [dE]- Parton energy loss, which 
can reproduce the flat, suppressed Raa in single particle spectra [ 3j, fails to reproduce 
the required V2 at intermediate px [ HI • The complicated interplay between soft and 
hard physics (for example, recombination, flow, and energy loss) at intermediate pt makes 
it critical to measure V2 at high pt, where we expect the only contribution to be from 
energy loss. We report here on the measurement of the 7r° azimuthal asymmetry, and 
discuss what it tells us about the source(s) of high pt f2- First we describe the method 
by which we measure Raa ^^s a function of angle with respect the reaction plane and V2 
for 7r°s in PHENIX. We then discuss two model calculations that help to shed light on 
the mechanisms giving rise to the high px V2- 



2. Measuring 7r° yields and t>2 in PHENIX 

During Run 4, PHENIX recorded 1.5 billion minimum bias Au+Au events at ^/s = 
200 GeV/c; the data presented here is taken from 1 billion of those events. For measuring 
photons and vr^s, we use the Electomagnetic Calorimeter (EmCal) The angle of 

the reaction plane is measured event-by-event using the Beam-Beam Counters. Pairs of 
clusters that pass 7 identification cuts are binned in angle with respect to the reaction 
plane {A(j) = — "^rp). A similarly binned mixed event background is then subtracted. 
The counts in the remaining peak are integrated in a ±2cr window, determined by a 
gaussian fit. Six bins in A0 are used from — 7r/2. 



*For the full list of PHENIX authors and acknowledgements, see Appendix 'Collaborations' of this volume 
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Figure 1. 7r° V2 and Raa'- |(a)| shows i?AA(PT) as a function of {Npart) for each A0 bin; 
the panels are for different bins. |(b)| the 7r° V2{pt), with each panel corresponding to 
a centrality bin. The grey bands indicate the systematic error due to the reaction plane 
resolution correction. 



To measure V2, we fit the A0 distribution as 1 + 2^2""' cos(2A0). The resulting V2 
parameter needs to be corrected for the reaction plane resolution, hence the designation 
^raw rjnj^g resolution a is determined for each centrality bin, and leads to the corrected 
value ^2°^^ = f^'^/cr. The yields as a function of A0 can then be corrected with a factor 

, l + 2<^^cos2A0 

f = -■ (1) 

1 + 2t;5"^cos2A0 ^ ^ 

3. Results and Discussion 

To obtain Raa{^4>), we exploit the fact that the ratio of of the yield at a given A0 to 
the inclusive yield is equivalent to the ratio of the angle-dependent Raa to the inclusive 
Raa- Thus mulitplying these relative yields by an inclusive measured Raa, we have: 

i?AA(A0) = Yield(A0) /Yield x Raa (2) 

The RAAi^(l>,PT) as a function of {Npart) is shown in Figure [T(a)j It is clear that there 
is non-trivial substructure to the angular dependence of the Raa, and that it varies with 
centralities. This feature is emphasized by plotting the data on a semi-log scale, showing 
that the Raa behaves differently in different A0 bins. 

The resulting vr^ f 2 is shown in Figure |l(b)| For the first time we observe f 2 up to 10 
GeV/c, and we see a clear decrease at high pt, then a leveling off to a finite value. 

To gain insight into the V2 mechanisms at work at high p^, we turn to models. We 
compare the f 2 to two models, an Arnold- Moore- Yaffe (AMY) calcuation [ [7] done 
by Turbide et al. and the Molnar Parton Cascade (MFC) model [|H]. Figure |21 shows 
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Figure 2. Comparison of 7r° f2 with models. The top three panels show the AMY calcu- 
lation with data for three centralities. The bottom two panels compare two centralities 
with the 6 = 8 fm calculation of the MPC model. 



calculations from these models, plotted alongside data for similar centralities. The AMY 
calculation contains energy loss mechanisms only, and we see that the data appear to 
decrease to a value at high px that is consistent with this model; the agreement is most 
striking in the 20-30% bin. 

The MPC model has a number of mechanisms in it, including corona effects, energy 
loss, and the ability to boost lower px partons to higher px (a unique feature). The 
calculation shown in Figure El does a better job of reproducing the overall shape of the 
V2, though it is systematically low. It is important to note that this calculation is done 
for one set of parameters, so it should be very interesting to see if the MPC can better 
reproduce the data for a different set of parameters. 

The prevailing mechanism for high pt f 2 is parton energy loss in the dense medium. If 
this is true, the Raa should be sensitive only to the geometry of the collision. To test this 
behavior, we seek to reparamaterize the two handles we have on geometry (centrality, or 
collision overlap, and angle of emission) into a single parameter, an effective path length 
which we will refer to as Lg^f . Details of the calculation are described in In essence, it 
is propotional to the parton-density weighted average of the length from origin to edge of 
an ellipse. We also perform a Glauber Monte Carlo sampling of starting points to account 
for fluctuations in the location of the hard-scattering origin of the particles' paths. 

The result of plotting Raa for all centralities and angles vs. L^jj is shown in FigureEl If 
the observed Raa arose from only geometric effects, we would expect the data to exhibit 
a universal dependence on Lgff. For low px, this is clearly not the case; something more 
than just energy loss is taking place there. However, when the pt reaches 7 GeV/c and 
above, the Raa data do indeed appear to have a dependence on a single L^ff curve. 
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Figure 3. Raa{^4>^Pt) vs. L^jj. The panels correspond to different pt ranges. The colors 
represent the centrality bins: magenta is most peripheral while black is most central. 



4. Summary 

We have presented the first measurement of high px f 2 for 7r° and charged hadrons. 
It is now clear that the V2 at pt > 7 GeV/c decreases to a small but non-trivial value. 
Comparison with models suggest that the dominant mechanism at work at high pt is 
energy loss. In addition, we have also presented the first measurement of 7r° Raa as a 
function of angle with respect to the reaction plane. The Raa{^4>,Pt) exhibits interesting 
angular substructure. Furthermore, when the Raa data are plotted as a function of an 
effective path length through the medium, they exhibit a universal behavior at px > 
7 GeV/c. 
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